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ABSTRACT 

In this paper we derive the Spitzer 7i?^Cband period-luminosity (P-L) relations for 
the Small Magellanic Cloud (SMC) Cepheids, by matching the Spitzer archival SAGE- 
SMC data with the OGLE-III SMC Cepheids. We find that the 3.6/im and 4.5/im band 
P-L relations can be better described using two P-L relations with a break period at 
log(P) = 0.4: this is consistent with similar results at optical wavelengths for SMC P-L 
relations. The 5.8/im and 8.0/im band P-L relations do not extend to sufficiently short 
periods to enable a similar detection of a slope change at log(P) = 0.4. The slopes of the 
SMC P-L relations, for log(P) > 0.4, are consistent with their LMC counterparts that 
were derived from a similar dataset. They are also in agreement with those obtained 
from a small sample of Galactic Cepheids with parallax measurements. 

Subject headings: stars: variables: Cepheids — distance scale 
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1. Introduction 

A precise and accurate measurement of the Hubble constant at the ^ 2% l evels is important 



for modern precision cosmology (for example, see iTegmark et al 



Jackson|[2007l : loihnglboOTl : lGreenhili[2009l : IXomatsu et ahlboiol : 



2004 



Hu 



Larson et al. 



2005:lMacri et al 



2010 



2006 



Mortonson et al. 



2010l ). This will be feasible after the launch of the James Web b Space Telescope fJWST), which will 
be ca pable of routinely observing Cepheids in distant galaxies ( Riess et al.ll2009l : lFreedman Sz Madore 
2010l ) at near- and mid-infrared wavelengths. In preparation for future distance scale work, it is 
necessary to derive the Cepheid Period-Luminosity relations (hereafter P-L relations, also known as 
the Leavitt Law) at mid-infrared (~ 3^m to ~ 8/im) wavelengths. Some adv antages and potentia l 
problems of the mid-infrared P-L relations were summarized, for example, in lNgeow et al.l (J2009al . 
and reference therein) . For example, the contribution of extinction to the error budget of the Hubble 



constant becomes negligible in these bands . Motivated by th i s, the Spitze r IRA C hand (h e reafte r 
IRA C band) P-L relations were derived in iNgeow &: Kanbuii (J2008al ) and iFreedman et al.1 (|2008l ) , 
by matching the known Cepheids in the Large Magel lanic Cloud (LMC) to the archival SAGE 
dataset (Surveying the Agents of a Galaxy's E v olutionjMeixner et al.ll2006). Th e IRAC hand P-L 
relations were further refined in iNgeow et al.1 (J2009bl ) and iMadore et al.l (|2009l ). respectively, by 
using the SAGE Epoch 1 and 2 data. Besides the LMC Cepheids, the IRAC hand P-L relations 
were also derive d from a sample of G alactic Cepheids that possess independent distance estimates 
from literature ([Marengo et al.ll2010l ). 



The main purpose of this Paper is to extend the IRAC band P-L relations to the Small 
Magellanic Cloud (SMC), which has a much lower metallicity than both the LMC and our Galaxy. 
The slopes of the P-L relations in the IRAC bands are expected to be insensitive to metallicity 
(jFreedman et al.l l2008l ) . An empirical determination of the IRA C band P-L relation derived from 
SMC Cepheids will provide a critical test of this assumption. 



The Data 



Recently, the Optical Gravitational Lensin g Experiment has relea sed a catalog of SMC Cepheids 
from its third phase of operation (OGLE-IIL ISoszvhski et al.ll20ld ). This contains 2626 funda- 
mental mode Cepheids. To derive the SMC ffl^Cband P-L relations, we matched the OGLE- 
III SMC Cepheids to the pubhcly released Epoch 1 SAGE-SMC (SAGE-SMC: Surveying the 
Agents of Galaxy Evolution in the Tidally-Disrup ted, Low-Metallicity Small Magellanic Cloud, 
Gordon fc SAGE-SMC Spitzer Lepcv TeaiaWmm IRAC hand data ( version 518.0.2). To be 
consistent with previous studies (JNgeow &: Kanbuii l2008al : iNgeow et al.l l2009bl ) , we only use the 
SAGE-SMC Archive data (hereafter SAGE-SMC data), which contains ^ 1.28 millions sources. 

We matched the SAGE-SMC sources to the OGLE-III SMC Cepheids by using a search radius 
of 2". The number of matched SAGE-SMC sources and the mean separation (A) from the input 
OGLE-III SMC Cepheids in each 7i?^C bands are summarized in Table [TJ From this table, it 
can be seen that more than 95% of the matched SAGE-SMC sources are located within 1", or 
within one pixel in the ZR^Cband images (the pixel scale for /Rj4 C instrument is ~ 1.2"/pixel, see 
Fazio et al.ll2004l ). from the OGLE- III SMC Cepheids. Extinction is ignored in this Paper, because 
it is negligible in the /i?^C bands (JFreedman et al.ll2008l ; iNgeow et al.ll2009bl ). 



3. The SMC P-L Relations in the IR^C Bands 



Figure [J displays the P-L relation in the four IRAC hands using all of the matched SAGE- 
SMC sources. This Figure presents a clear correlation between the IRAC hand magnitudes and the 
Cepheid pulsation period. To remove the outliers shown in the Figure, we employed an iterative 
2.5o" rejection algorithm when fitting the P-L relations to the data. This is the same procedure used 



Table 1. Number and mean separation of the matched SAGE-SMC sources. 



Band 



A^„ 



< A>=' 



Std. Dev.'' Fraction within 1" 



3.6/im 


2567 


0.263 


0.214 


98.75% 


4.5/im 


2545 


0.264 


0.216 


98.70% 


5.S^J,m 


699 


0.265 


0.271 


97.00% 


S-O^m 


404 


0.284 


0.290 


96.29% 



^A is the separation, in arcsecond, between the matched SAGE- 
SMC sources and the OGLE-III SMC Cepheids. 

''The standard deviation of the mean. 

'^The fraction of matched SAGE-SMC sources within 1" radius 
from the OGLE-III SMC Cepheids. 
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Fig. 1. — The IRAC band P-L relations for all the matched SAGE-SMC sources to 
the OGLE-III SMC Cepheid, using a search radius of 2". The data were separated in 
two groups: one group with separations (A) between the Cepheid and matched sources 
that are greater than 1.2" (open circles), and another group with A smaller than 1.2" 
(small dots). The dashed horizontal lines are the faint limits in the respected bands, 



adopted from the SAGE-SMC document (see Ihttp : //data, spitzer . caltech . edu/popular/ 



sage-siiic/200911ir_enhanced/documents/sage-smc_delivery _nov09.pdf). P is the pulsating 
period for Cepheids in days. Error bars are omitted for clarity. 
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Fig. 2. — Two examples of the abnormally bright outliers shown in Figured! The left panels are the 
postage-stamp imag es for the two OGLE -III bands and four ZiZ^C bands. The VI band images 
were adopted from lUdalski et al.l (J2008l ). and the IRAC images were downloaded from Spitzer 
Science Center data archive tool, Leopard. The circles highlight the location of the Cepheid, which 
have a radius of 2" . The large filled circles in right panels show the location of this Cepheid in the 
respected P-L relations. These two examples do not have any nearby neighboring sources within 
1.2" from the Cepheid. 



OGLEIII-SMC-CEP-1251: 

Neighbouring star: I = 18.89, no V band 
pinotometry and (V - I) color 
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Neighbouring star: I 
(V-l) = 0.59 
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Fig. 3. — Same as Figure [21 but for two Cepheids that have a nearby neighboring source within 
1.2" from the Cepheid. The location of the neighboring source is indicated with smaher circles 
(with radius of 1") in the left panels. 



in iNgeow &: Kanbuij (|2008al ) and iNgeow et al.1 (|2009bl ) . The abnormally bright outliers shown in 
Figured] could be blended with a nearby red companion. To find out if there are any nearby sources 
around these outliers, we searched for neighboring so urces within 1 . 2^^ ra dius from the Cepheids' 
location, using the OGLE-III SMC photometry map (jUdalski et al.ll2008l . which has a pixel scale 
of 0.26"/pixel). Only about 20% of of these outliers were found to have one faint neighboring 
source around the Cepheids (and few of them have more than one neighboring sources). Figure 
[2] shows two examples of the outliers that do not have any neighboring source; while Figure [3] 
presents two examples that display a neighboring source. These abnormally bright outliers could 
be due to the presence of circumstellar dust envelopes that would cause the observed mid-infrared 
excess. The circumstellar e nvelope could be formed from Cepheid mass-loss activity (see, for 
example. iNeilson et al.ll20ld . and reference therein). A detailed investigation of these outliers and 
their mass-loss activity is beyond the scope of this Paper, but will be addressed in future papers. 
Nevertheless, it is clear that these outliers should be removed from the sample. 

Figure [H also suggests that a period cut is needed for the 5.8/xm and S.O^um band data, as the 
data approaches the faint magnitude limit for the short period Cepheids. We therefore applied a 
range of period cuts to the data. Figure H] shows the slopes of the P-L relations for all four IRAC 
bands as a function of period cut. We discuss the implications of the period cuts and the resulting 
P-L relations in the following sub-sections. 



3.1. The 3.6/im & 4.5/im Band P-L Relations 

Figure [T] implies that period a cut may not be needed for the 3.6/im and 4.5/xm band data. If a 
period cut is not applied to the data, the resulting P-L relations are: ms.g/^m = — 3.372(ib0.012) log(P)+ 
16.570(±0.006) {N = 2378) and m^,^^^ = -3.338(±0.014) log(P) + 16.508(±0.007) {N = 2389). 
However, Figure U] points to the fact that the slopes of the P-L relations in these two IRAC hands 
gradually changes from a steep slope (~ —3.35) to a more shallow slope (~ —3.2) when the the 
adopted period cut, log(Pcut)) is less than ~ 0.4. In the optical b ands, a change of slope for SMC 
P-L relations at 1ok(P) ~ 0.4 has been reported in the literature (IBauer et al.lll999l : lUdalski et al. 



19991 : ISharpee et al.ll2002l : ISandaee et al.ll2009l: ISoszvhski et al. 



20ld'). This slope change is postu- 



lated to occur due to evolutionary effects (jBecker et al.lll977l : iBaraffe et al.lll998l ). Therefore, we 
separated out the 3.6/im and 4.5//m band data at log(P) = 0.4. For Cepheids with log(P) < 0.4, 
the resulting P-L relations are: 



ms.e^m = -3.546(±0.049)log(P) + 16.617(±0.011), 

^3.6Mm = 0.183, iV3.6/.m = 1478, 

m4.5/.m = -3.414(±0.058)log(P) + 16.537(±0.013), 

fT4.5Mm = 0.214, iV4.5Mm = 1481. 



While for Cepheids with log(P) > 0.4, the resulting P-L relations are: 



ms.e^m = -3.225(±0.020)log(P) + 16.448(±0.015), 

o-s.e/xm = 0.163, iVs.e^m = 906, 
m4.5f,ni = -3.177(±0.021)log(P) + 16.371(±0.016), 

0-4.5Mm = 0.172, iV4.5Mm = 906. 



The difference between the long (log[P] > 0.4) and short (log[P] < 0.4) period slopes are at the 
~ 60" and ~ 4cr level for the 3.6//ni and 4.5/ini band P-L relations, respectively. Figure [5] presents 
the residuals of the fitted P-L relations as a function of period when the long period P-L relations 
are used to fit all data. The residuals show a clear trend at the short period end. To furt her verify 
that there is a change of slope at log(P) = 0.4, we applied the F-test (for example, see IWeisberg 
1980) to the combined data of the long and short period Cepheids (after removing the outliers). 
The F-test results show that the 3.6;uni and 4.5/xm band P-L relations are non-linear at the break 
period of log(i-') = 0.4, with Fsg/^m = 39.2 and -F4.5/xm = 33.1, respectivelju. 

For Cepheids with log(P) > 0.4, we also tested the non-linear P-L relations at a break period 
of log(. P) = 1.0 (which have been found to be non-linear for LMC P - L relations in the opti 
cal, 



see 



Kaiibur fc Ngeow l2004l : ISandage et al. 



2004 



Ngeow et al. 



2005 



Kanbur fc Ngeow 



200e 



Ngeow et al.ll2008bl . l2009bl . and reference therein). The F-test results imply that the P-L relations 
are linear (with -Fs.fi j/m = 0.35 and F4^5/;,m = 0.07). This is consist ent with tests on the LMC IRAC 



band P-L relations (jNgeow k Kanbuill2008al : iNgeow et al.ll2009bl ). 



3.2. The 5.8/im & 8.0/xm Band P-L Relations 

It is clear from Figure [T] that the shallow slopes for 5.8;uni and 8.0^ni band P-L relations 
with log(Pcut) < 0.5, as shown in Figure HI are due to the lack of Cepheids at the short period 
end, where the magnitudes approach the faint limit in these two IRAC bands. A period cut of 
log(Pcut) > 0.5 is necessary to remove the bias due to the incomplete data at the faint magnitude 
end. Inspecting Figure S] suggests that the slope of the 5.8/im and 8.0^m band P-L relations begin 
to stabilize at log(Pcut) ~ 0.7 and log(Pcut) ~ 0.8, respectively. Hence, we have adopted a period 
at log(Pcut) = 0.7 and log(Pcut) = 0.8 when deriving the P-L relations in these two bands. The 
resulting P-L relations are: 



"l5.8^m 



-3.227(±0.037)log(P) + 16.381(±0.038), 



^Recall that for a large number of data points, F ~ 3 at 95% confident level. Therefore, _F > 3 indicates that the 
null hypothesis of a single-line regression can be rejected. 
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Fig. 4. — The slopes of the P-L relation as a function of adopted period cut (log[Pcut])) where only 
the data with log(P) > log(Pcut) were used in deriving the P-L relations. 
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(T5.8^lm = 0.180, iVs.S/im = 334, 

ms.Ofim = -3.253(±0.048)log(P) + 16.400(±0.053), 
O'S.Ofj.m = 0.189, Ns.O^m = 227. 



The ratios of the number of Cepheids with log(P) > 1.0 and log(P) < 1.0 are 1 : 2 and 1 : 1 in 
the 5.8//m and 8.0//m band, respectively, after the period cut. Therefore the F-test can be apphed to 
the data in order to test for non-hnearity at log(P) = 1.0. The F-test resuhs show that both of the 
5.8//m and 8.0/xm band P-L relations are linear, wit h Fq.8/;.m = 0.06 and -Fs.n/Arp = 0.08 respectively . 



This result is also consistent with the LMC results (JNgeow fc Kanbuiil2008al : iNgeow et al.ll2009bl ) 



4. Comparison of the P-L Slopes 

In Table [21 we compare the ZR^Cband P-L slopes currently available for our Galaxy, LMC 
and SMC. We aggregate the slopes into two groups: one with shallow slopes and another with 
steep slopes. In gener al, slopes within a gr oup are consistent with each other. In contrast, the 
LMC slopes derived in iNgeow et al.l (J2009bl . the LMC2) and the SMC slopes disagree with those 
in the "steep slope" group at the 3a level or more. If the IRAC band P-L slopes were indeed 
shallow, then the slopes from all three galaxies would be consistent (despite the large uncertainties 
for GAL3 slopes shown in Table [2]). This suggests that the metallicity will not strongly affect the 
IRAC hand P-L slopes. On the other hand, if the Galactic and LMC P-L slopes were steeper (for 
those listed as GALl/2 and LMCl in Table [2]), our SMC slopes will challenge the assumption that 
these slopes should be insensitive to metallicity. 

Does crowding affect the derivation of the SMC P-L slopes presented in this Paper? We tested 
this by dividing our SMC Cepheids into two groups: one group for Cepheids located in the central 
"bar" region (with higher stellar density), and another group for Cepheids located on the outskirts 
of the "bar" region. The derived P-L relations from these two groups are consistent with each other. 
We performed another test by searching the nearby companion sources within 1.2" of the respected 
Cepheids, using the OGLE-III SMC photometry map mentioned previously. The Cepheid samples 
were grouped into two groups that either did or did not contain the nearby sources. Again, the 
slopes derived from these two groups are in good agreement. We therefore believe crowding does 
not seriously affect our results. 

The discrepancy between the two L MC slopes given in Table [2] may be due, in part, to the 



different samples used in b oth studies. iMadore et al.l (l2009l ) derived the P-L relations using 70 
LMC Cepheids drawn from lPersson et al.l ( 20041 ) . Most of the i n had p eriods longer than 10 days. 



In contrast, the LMC P-L relations presented in INgeow et al.l (J2009bl ) were derived based on the 



1 pt 

m 



~ 1800 OGL E-III LMC Cepheids (ISoszvhski et al.l 2008 ). which have a distribution peaking at 
log(P) ~ 0.5 (|Soszvhski et ahlboioluNeilson et alT(|201ol ) showed that the LMC I RAC hand P-L 



slopes, based on the OGLE-III Cepheids, became comparable to the slopes given in lMadore et al 
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(J2009l ) when only the long period (log[P] > 1) Cepheids were used. 



For the Galactic ffl^Cband P-L relations, the P-L slopes derived from the infrared surface 
brightness (IRSB) technique (GALl and GAL2 in Table [2|) were steeper than the slopes derived 
from Cepheids that have parallax measurements (GALS in Table [2]). The large uncertainties of 
the GAL3 P-L slopes (due to t he combination of a small number of data points in the sample 
and less precise photometry, see iMarengo et al.1 120101 ) cause all three sets of Galactic P-L slopes 
to be consistent w i th ea ch other. It is worthwhile to note that the theoretical P-L slopes derived 



m iMarengo et al.l (|20ld ) agree well with the GALS P-L slopes. If the uncertainties of the P-L 
slopes derived from the parallax distances are the same as those from the I RSB techniques, then 
the slopes between them will disagree at the ^ 3.5a level. IMarengo et al.l (|20ld ) suggested that 
the projection factor-period relation used in IRSB technique may need further refinement. 

Finally, we have derived the AKARI N3 hand (at ~ 3 / /m) P -L relation for the LMC Cepheids, 
which is independent of the SAGE data, in iNgeow et al.l (|20ld ). The slope of the A^S band P-L 
relation is — S.25 ± 0.05. This agrees well with the 3.6/xm band P-L slopes in Table [2] for the slopes 
in the "Shallow Slope" group. 



Conclusion 



In this Paper we have derived SMC P-L relations in the IRAC hands, by matching the archival 
SAGE-SMC data to the latest SMC Cepheid catalog from OGLE-III. We have found, for the first 
time, that the SMC P-L relations show a change of slope at log(P) = 0.4 in the 3.6/um and 4.5^m 
band, similar to their optical counterparts. Due to the incompleteness at the short period end, such 
a change of slope cannot be confirmed to exist in the 5.8//m and S.O^m bands. Future observations 
with JWST may help to determine the change of slope in these two bands. 

The slopes of the SMC IRAC band P-L relations were fou nd to be around — S.2, which is 
consistent with the LMC slopes found in iNgeow et al.l (l2009bl') and th e Galactic P-L relations 
derived from Cepheids with parallax measurements ([Marengo et al.ll201G ). This indicates that the 
slopes of the P-L relations in the IRAC bands are insensitive to metallicity. Howe ver, the SMC 



slopes disagree with the steeper slopes for the LMC and Galactic counterparts from lMadore et al. 



(J2009l ) and those from the IRSB technique, respectively. Future observations such as the parallax 
measurements from Gaia for Galactic Cepheids and/or JWST observations of Magellanic Cloud 
Cepheids may help to resolve this discrepancy. 
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Table 2. Comparison of the P-L slopes in the /i?^C bands. 



Galaxy 



3.6/xirL 



4.5/im 



5.8/im 



8.0/im 



Reference^ 



Note 



Steep Slope 



GALl 


-3.54 ±0.04 


-3.45 ±0.04 


-3.51 ±0.03 


-3.57 ±0.03 


(1) 


GAL2 


-3.47 ±0.06 


-3.38 ±0.06 


-3.44 ±0.06 


-3.46 ±0.06 


(1) 


LMCl 


-3.40 ±0.02 


-3.35 ±0.02 


-3.44 ±0.03 


-3.49 ±0.03 


(2) 



29 GAL Cepheids with "old' 
28 GAL Cepheids with "new 
70 LMC Cepheids 



IRSB distances'' 
" IRSB distances'" 



Shallow Slope 



GAL3 


-3.16 ±0.22 


LMC2 


-3.25 ±0.01 


SMC 


-3.23 ±0.02 



-3.06 ±0.23 -3.10 ±0.23 -3.16 ±0.18 (1) 

-3.21 ±0.01 -3.18 ±0.02 -3.20 ± 0.04 (3) 

-3.18 ±0.02 -3.23 ±0.04 -3.25 ± 0.05 (4) 



8 GAL Cepheids with parallax distances'' 
~ 400 to ~ 1630 LMC Cepheids 
~ 220 to ~ 900 SMC Cepheids 



''Reference: qi lMarengo et all jioT^); (2') lMadore et alj||2009h : (3') lNgeow et al.l ll2009bl 'l: (4) This Paper for Cepheids with log(P) > 0.4. 

''"Old" and " new" distances referred to the old and new projection factor-period relations adopted in IRSB technique, respectively. See 
iMarengo et al.l 1120101) for more details. 



